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We focus on the probability distribution function (pdf) P( Ay; 7) where A7 are the measured 
strain intervals between plastic events in an athermal strained amorphous solids, and 7 measures 
the accumulated strain. The tail of this distribution as A7 —>• 0 (in the thermodynamic limit) scales 
like A7’ 7 . The exponent 77 is related via scaling relations to the tail of the pdf of the eigenvalues of 
the plastic modes of the Hessian matrix P(A) which scales like X e , 7 = (6 — l)/2. The numerical 
values of 7 or 0 can be determined easily in the unstrained material and in the yielded state of plastic 
flow. Special care is called for in the determination of these exponents between these states as 7 
increases. Determining the 7 dependence of the pdf P( A7; 7) can shed important light on plasticity 
and yield. We conclude that the pdf’s of both A7 and A are not continuous functions of 7. In slowly 
quenched amorphous solids they undergo two discontinuous transitions, first at 7 = 0 + and then at 
the yield point 7 = 7 Y to plastic flow. In quickly quenched amorphous solids the second transition 
is smeared out due to the non existing stress peak before yield. The nature of these transitions and 
scaling relations with the system size dependence of (A7) are discussed. 


I. INTRODUCTION 

The response to external strain in amorphous solids 
as the strain 7 is increased in an athermal quasistatic 
manner (AQS) does not result in a smooth differentiable 
stress-strain relationship. Rather, it results in a sequence 
of reversible strain increases A71, A72, A73 • • • followed 
in each case by an irreversible plastic event in which a 
local, sub-extensive or extensive group of particles dissi¬ 
pate energy and the system then falls into a new inherent 
state. After each such event A7, the stress has changed 
by A Oi. An ensemble of sufficiently many independent 
realizations of this process at any value of 7 will then give 
rise to a probability distribution function (pdf) P(Ay; 7). 
The properties of this pdf, and in particular of the scaling 
form of the tail of this pdf as A7 — > 0 were the subject 
of focused recent research. The tail of the distribution 
is relevant for the thermodynamic limit since the inter¬ 
vals A7i shrink when the system size increases as shown 
below, Eq. (Ell¬ 
in a pioneering paper in 2010 lj it was shown that for 
7 = 0, i.e. after a quench from the liquid state to the 
amorphous state at T = 0 the tail of the pdf appears to 
scale like 

lim P(A 7 ,7 = 0) ~ (Ay)* 7 . (1) 

A7—>-0 

The exponent 77 was measured there to be 7 ft 0.6 in 2- 
dimensions without any claim for universality; in partic¬ 
ular this exponent may depend on the quench rate from 
the liquid to the solid. Besides the presence of a scaling 
tail it was shown that the whole pdf agrees extremely 
well with the Weibull distribution, giving rise to an in¬ 
teresting scaling relation to the system size dependence 
of the mean value of A7 at 7 = 0, 

(A 7 ) ~ N& , /3 < 0 . (2) 

The scaling relation follows from the theory of extreme 


value statistics, and with the Weibull form it reads [lj 

, = “( 1+ £) ■ (3) 

Indeed, the exponent /3 was found to be in perfect con¬ 
sistency with 77 Ri 0.6, i.e. /? ft —0.62. The method of 
measurement employed many independent realizations of 
quenched systems at 7 = 0. It is important to stress that 
for large systems the mean value (A7) tends to decrease 
with system size according to Eq. ©■ The scaling pre¬ 
sented in Eq. © refers always to values of A7 -C (A7). 

Other statistically stable scaling relations can be ob¬ 
tained in the steady state plastic flow after the plastic 
yield. Here we have many repeated events in stationary 
conditions, allowing us to measure carefully the statis¬ 
tics of energy drops A Ui, stress drops Acq and strain 
intervals A7*. It was found numerically in [2] and then 
analytically in [l| that 

(A U) = eN a , (A a) = sN 0 , (A 7 ) = , (4) 

with a = 1/3 and /3 = —2/3 as exact universal results. 
On the other hand it was shown that the relation © is 
no longer valid in the plastic flow state. In fact there 
77 = 0 and the pdf has no relation to the Weibull distri¬ 
bution. This major difference was ascribed to the exis¬ 
tence of sub-extensive plastic events in the plastic flow 
steady-state. These sub-extensive events are correlated, 
destroying the statistical independence of repeated real¬ 
izations in the quenched 7 = 0 state which is necessary 
for the scaling laws ©. 

The changes in the statistical properties as they are 
reflected in the pdf of the strain intervals A7 raises an 
interesting challenge. Can one follow these changes to 
understand plasticity in amorphous solids, and in partic¬ 
ular to study the yielding transition to the steady state 
plastic flow. This challenge was picked up in particular 
in Refs. 1,1- The upshot of these studies is that there 
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exists a continuous change, as a function of the external 
strain 7, of the values of the scaling exponent, from 7 = 0 
to the plastic flow regime. If correct this would supply 
an extremely rare occasion of a continuous change in a 
scaling exponent as a function of a control parameter Q. 
The aim of this paper is to study carefully this interest¬ 
ing proposition. We end up offering a different point of 
view in which we stress two discontinuous transition, one 
at 7 = 0 + and the other at the yield point 7 = 7 Y - It 
should be pointed out however that in this paper we dis¬ 
cuss only the measurable strain intervals between plastic 
events and the lowest eigenvalues of the Hessian matrix 
that vanish during the approach to a plastic events. In 
Refs. @,0 there appear distributions of variables that 
are not readily measurable in strained amorphous solids, 
and a full comparison of predictions is not always easy. 
Nevertheless our conclusion that the yield is a discontin¬ 
uous transition (at least in slowly quenched amorphous 
solids) is in accordance with recent studies that have used 
very different methods of investigation j@]. 

In Sect. El we present some numerical estimates of the 
exponent /? and of the pdf P( A7; 7) for a range of values 
of the external strain 7. To interpret properly the nu¬ 
merical results we turn in Sect. nn to theoretical consid¬ 
eration. These considerations are based on relating the 
scaling properties of P{ Ay;7) to the scaling properties 
of another pdf, P{ A; 7), which describes the distribution 
of eigenvalues of the plastic modes of the Hessian matrix, 
precisely those eigenvalues that approach zero before a 
plastic event is taking place. The separation into De¬ 
bye modes and plastic modes was justified in 00 , and 
see also 0- In Sect. ED we present a Fokker-Planck 
equation for the development of P{ A; 7) as a function of 
increasing 7. While we cannot vouch for the exactness 
of this equation, we argue that the scaling properties of 
this equation are robust and can be trusted to provide 
the correct scaling exponents for the tails of P( A; 7) and 
P(A75 7). In Sect. HVl we offer a summary and conclu¬ 
sions. 


II. DIRECT NUMERICAL RESULTS 

In this section we describe direct numerical measure¬ 
ments of the exponent /3 and 7. To generate data we have 
employed a 50-50 binary mixture of point particles inter¬ 
acting via Lennard-Jones potential in 2 -dinrensions. The 
parameters of the model can be found in Ref. a The 
system is quenched from a high temperature liquid to a 
target temperature T = 0.001 in Lennard-Jones units. 
We use two quench rates. In the first set of simulations 
the quench is “infinitely fast” in the sense that we use 
conjugate gradient energy minimization. In the second 
we reduce the temperature at a rate of 10 ~ 5 in Lennard- 
Jones units (cf. Ref. 0 ), and we refer to it below as the 
“slow” quench. The mechanical properties of the differ¬ 
ently quenched systems differ. The slow quench shows 
a distinct stress peak when 7 is increased, followed by 


a yield toward the steady state plastic flow. The fast 
quench results in a gradual increase in stress towards the 
plastic flow steady state without a stress peak. 

After preparing the system in athermal condition we 
follow the standard AQS protocol to strain the system 
in simple shear strain. We carefully measured the strain 
intervals A7* that occurred between plastic events. We 
back-tracked our simulations to increase the precision of 
this measurements (cf. Appendix of 0 )- Every such sim¬ 
ulation was repeated using 1000-1500 freshly quenched 
independent realizations. We have measured the distri¬ 
bution of A74 and the mean value (A7) as a function of 
N for different system sizes. While there is no ambiguity 
about the latter quantity at 7 = 0, for higher values of 7 
we need to collect data within a bin of values of 7. We 
opted to do so in bins of size 0 . 01 . Thus where we report 
below a value of 7, say 7 = 0.02 one needs to interpret 
that as data collected from all events that occur between 
7 = 0.01 and 7 = 0.02. 

In Fig. [T] we show some representative results for the 
measurement of the exponent / 3 . In the upper panel the 
scaling law ([2]) is demonstrated for 7 = 0, whereas in the 
lower panels we show the result for 7 = 0.04 and 7 = 0 . 1 . 
The first one is in agreement with the previous estimate 
of ~ — 0 . 62 . The last one is beyond the plastic yield 
and is in agreement with the exact prediction /? = — 2 / 3 . 
The middle panel is an example of the new numbers ob¬ 
tained here for the first time for intermediate values of 
7. An over-all impression of the resulting picture can be 
obtained if we assume that the scaling law m is obeyed 
throughout the range of 7. We can then use the measure¬ 
ment of /3 to plot — (1 +jj) as a function of 7, see Fig. [ 2 ] If 
one believes that this represents the exponent 77 accord¬ 
ing to ©, then this continuous variation as a function 
of 7 is somewhat reminiscent of the dependence of this 
exponent on 7 as predicted in Refs. 0 , 0 | (This exponent 
is denoted as 9 there). In fact we will argue now that 

(i) the scaling law © is valid only as long as 7 < 7 Y . 

(ii) The measured values of /3 for 0 < 7 < 7 Y suffer from 
severe finite size and cross-over effects. In fact the anal¬ 
ysis presented below indicates that 77 = 0 and /3 = — 1 
for all the values of 7 in the range 0 < 7 < 7 Y , with a 
transition to /? = — 2/3 after the plastic yield transition. 
The transition is expected to be sharp in the case of the 
slow quench, but much more smeared out in the case of 
the fast quench. 

To see that the scaling law © fails for 7 > y Y it 
is sufficient to examine the pdf’s P(In Fig. [ 3 ] 
we present these distribution functions for the slowly 
quenched amorphous solid for a range of 7 values start¬ 
ing with 7 = 0 and ending with 7 = 0.09 which is 
already beyond the plastic yield which for the present 
slowly quenched system is estimated at y Y ~ 0 . 06 . It is 
clear that for 7 = 0 the pdf has a perfect Weibull form 
with r] ss 0.6 in agreement with [3 ss — 0 . 62 . But for all 
7 > 0 the pdf has a flat tail as A7 —> 0 , indicating a 
value of 77 = 0 . In the plastic flow state, at 7 > 7 Y , 
where we measure /3 = — 2 / 3 , the scaling law ([ 3 |) is not 


3 





FIG. 1 . The scaling law (j 2 j demonstrated for 7 = 0 (upper 
panel), 7 = 0.04 (middle panel) and 7 = 0.1 (lower panel). 
We propose that in the range 0 < 7 < 7 y the value of (1 in 
the thermodynamic limit is /3 = — 1 . 


obeyed. It was proposed in Ref. . 1 ] that this is due to 
the sub-extensive plastic events that occur in the flowing 
state. The arguments leading to the scaling law © fail 
under these conditions. 

The question is then whether the scaling law © holds 
in the range 7 < y Y . We will argue that the answer is in 
the affirmative. We read 77 = 0 from the results of Fig. © 
and we propose that the correct value of /? is /? = — 1, in 
agreement with the scaling law ©. This is equally valid 
for the slow and the fast quench (see the results of the 
pdf’s for the slow quench in Fig. ©. The results shown 
in Fig. [2] should be read with this in mind; due to system 
size and crossover effects the sharp transition at 7 = 0 + 
is smeared out, and so is the transition near 7 Y . Clearly, 
in the case of the fast quench the transitions are even 



FIG. 2 . A plot of — ( 1 T ^) obtained from the direct measure¬ 
ment of p as a function of 7. Results for the slow quench are 
shown in circles and for the fast quench in squares. We argue 
in this paper that these results suffer from severe finite size 
effects and in reality this figure should be replaced by Fig. [G] 


Slow Quench, N = 4000 




Ay 

FIG. 3 . Results for the slow quench. Upper panel: an exam¬ 
ple of the raw data for the probability distribution functions 
P( Ay;7) for 7 = 0 . 04 . Lower panel: the fits to the pdf’s 
P( A7;7) for a range of values of 7 from 7 = 0 to 7 = 0.9 
(beyond the plastic yield). Only at 7 = 0 the distribution is 
Weibull with 7 ~ 0 . 6 . For all other values of 7 we find 7 = 0 . 
This value is explained by the theory in Sect. m 
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FIG. 4. Results for the fast quench: the probability distribu¬ 
tion functions P(Ay; 7) for a range of values of 7 from 7 = 0 
to 7 = 0.9 (beyond the plastic yield). Only at 7 = 0 the dis¬ 
tribution is Weibull with 77 as 0.3. For all other values of 7 we 
find 77 = 0. This value is explained by the theory in Sect. m 


more smeared out. We propose that in the case of slow 
quench, in the first transition 77 tries to reach the value 
0 and in the second transition the scaling law © breaks 
down and /3 settles on the value - 2 / 3 . The theoretical 
predictions will be presented below in Fig. [G] after the 
considerations leading to these assertions become clear. 


III. THEORETICAL CONSIDERATIONS 


will allow the calculation of the tails of the distribution 
of the strain events P(Ayjy) at fixed 7 and the tail of 
the associated energy barrier distribution P £ (e; 7) at fixed 
7. We reiterate that the distribution P(A; 7) discussed 
here refers to the minimal eigenvalues that are exposed 
in the tail of the distribution as they tend to zero due 
to plastic processes. There may be other putative plastic 
modes that are not yet ready to express their potential 
instability and these are not taken into account here. A 
good example of these plastic modes are those discussed 
explicitly in Ref. Ed¬ 
it is very important to keep separate in one’s mind two 
sets of experiments. In the situation described above the 
strain 7 is fixed at some value and an ensemble of A77 
measured. A second experiment involves ramping up 7 
from zero and measuring consecutive reversible A7i fol¬ 
lowed by avalanching. The nth such event will occur at 
7 = £” =1 A 7,;. In such an experiment the stochastic dy¬ 
namics will prove to be very similar to particles drifting 
in a one dimensional flow together with added shot noise 
(with the eigenvalue A representing the particle ‘coordi¬ 
nate’ and 7 equivalent to time). In the thermodynamic 
limit, when the distribution of eigenvalues of the Hes¬ 
sian becomes a density, we will employ a Fokker-Planck 
equation for P(A; 7) which is similar in structure to a 
one dimensional flow for A with an adsorbing boundary 
at A = 0 due to the associated saddle-node instability 
followed by reinjection of the of a new eigenvalue with a 
distribution ji n ( A, 7) in order to conserve probability. Its 
generic form will be 

dP /8 7 = -d[v(X, 7)P] + DW&P/d A 2 + j in ( A, 7). (6) 


A. Relation to the Hessian matrix and its 
eigenvalues 

To understand the statistics of A7 we need to connect 
it to the statistics of the eigenvalues of the Hessian matrix 
El El- The latter is crucial for the understanding of 
atlrermal plasticity. Each plastic instability is due to a 
saddle-node bifurcation in which the lowest eigenvalue 

r\2 tt 

A i of the associated Hessian matrix Hij = q°.q t . of the 
potential energy U softens and becomes zero at some 
value of 7 = 7p. The modes whose eigenvalues become 
zero at some value of the external strain are referred to 
as plastic modes One should be aware that every 

disordered solid has also Debye modes whose eigenvalues 
remain roughly independent of the external strains. For 
the plastic modes, because we are dealing with a saddle- 
node instability, we know (and see below for details) that 
A7i ~ A? and also that the energy barrier £j scales as 
£j ~ A/. Analyzing the properties of an ensemble of the 
lowest eigenvalues of the plastic modes of the Hessian 
matrix at fixed strain 7 will yield (in the thermodynamic 
limit) a distribution P(A; 7). Knowing the tail of this pdf 
for A —> 0 , 


For understanding the scaling properties of P(A;7) it is 
crucial to derive the scaling form of v(A, 7) which we 
discuss next. 


B. The rate of change of the eigenvalues and the 
resulting scaling relations 


The central result that will be used in the present con¬ 
text is an “equation of motion” for the change in the 
minimal eigenvalues of the Hessian matrix upon increas¬ 
ing strain. This equation was derived in Ref. Q, and 
it turns out to be rather complex. Nevertheless, in the 
limit of A —> 0 this equation can be simplified to exhibit 
only the most singular contribution, correct up to less 
singular terms: 


dX a( 7 ) 

— ~---b less singular terms. 

<27 A 


( 7 ) 


Together with the boundary condition that A vanishes at 
some strain value 7 p, i.e. A(7 p) = 0 , we integrate Eq. © 
to obtain 


lim P(A; 7 ) 

A—>0 


( 5 ) 


A(q) = v / 2o(7p)(7p — 7) + higher order terms . (8) 
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We can also expand the projection of the potential energy 
in the direction that is becoming unstable for 7 —> "/ P 

U = U 0 + ±As 2 + lbs 3 + 0(s 4 ) . (9) 

This form implies that a saddle point exists at s* = — , 

of magnitude 


A E = U(s+) -U 0 ~ 


2A 3 

3 P 


4-^/2 

HT 


' ¥ (7P-7) 5 , (10) 


where we have used the solution Eq. [8] for the 7 depen¬ 
dence. Eq. © and Eq. m are the basic relations from 
which the scaling relations by which the strain increases 
and energy barriers are related to A 

A7 ~ A 2 (11) 

e ~ A 3 . (12) 

can be found. In particular Eq. m leads to the obvious 
scaling relation 


77 = (0-l)/2. (13) 

We expect the scaling relation m to be valid for all 
values of 7 and for all amorphous solids since it stems 
from the saddle-node nature of the plastic events which 
is generic. 


C. Consequences for the Fokker-Planck Equation 

The upshot of the preceding discussion is that we can 
estimate the most singular form of the ‘speed’ v(A; 7) in 
Eq. © as 

u(A; 7) « —a(7)/A + less singular terms . (14) 

We therefore rewrite Eq. © in the form 

dP/d 7 « o( 7 )d[P/A]/0A + D( 7 )d 2 P/dX 2 + j in (A >7 ) 

(15) 

with a probability flux 

J(A, 7 ) = -o( 7 )P/A - D( 7 )dP/d\ . (16) 

The diffusive term in the Fokker-Planck equation results 
from the “less singular terms” in Eq. m■ Note that we 
should seek solutions of this equation with an adsorbing 
boundary condition at A = 0, leading to 

P(A = 0,7) = 0, (17) 

but with a finite flux for any value of 7. We stress that 
we have added a non-singular (in A) diffusion term to the 
drift term in the Fokker Planck equation. The reason is 
that the theory leading to Eq. m does not indicate 
the existence of any other singular term in the “noise” 
on top of the leading, explicit 1/A term in Eq. (fl4l) . The 
most singular contribution is already there. For the same 



— isotropic 

— steady state 


FIG. 5. Plots of the isotropic Weibull distribution and the 
Steady State distribution assuming that the new eigenvalues 
are equally distributed between A = 0 and A = A max- The 
distributions have been chosen so that (A)i so = (X) 3S . 


reason also the injection term ji n { A,7) is taken as non¬ 
singular. 

Accepting the fact that for any value of 7 > 0 the 
flux must be finite, a direct consequence of the form of 
Eq. (gSD is that P(A;y) must start linearly in A in order 
to have a finite flux which is neither zero nor infinity: 

lim P( A; 7) ~ A . (18) 

A-s-0 

Using the scaling relation (1131) we conclude that 
77 = 0 for any value of 7 > 0. This result rational¬ 
izes entirely the data shown in Fig. [3j 

The same conclusion can be obtained directly by ana¬ 
lyzing the pdf in the steady state plastic flow. There both 
D("f) and ji n (X,"f) must become 7 independent, and we 
will assume a constant reinjection rate jm(A) = ji n be¬ 
tween 0 < A < X max in line with the avalanching process 
that occurs in the dynamics, while ji n {X) = 0 otherwise. 
We need to solve 


0 = a ss d[P/X]/dX + D ss d 2 P/d A 2 + j in ( A). (19) 

together with Eq. The solution is given by 

Pss(X) = y——[(A/A max ) ^ (X/Xmax) 2 ] (20) 
Amax 


as a quick substitution will confirm. It is properly nor¬ 
malized and (A) = X max /2. Choosing the same value for 
(A) for both the isotropic and steady state distributions 
yields Fig. [5] In addition, for this solution the probabil¬ 
ity flux J at the adsorbing boundary and the required 
reinjection rate ji n are respectively 


J = 


Jin — 


6(o> ss + 

Y2 

/x max 

6 {a ss + 2 D ss ) 

A 3 

m.n/r. 


( 21 ) 


The important thing to notice is that the conclusion 
Eq. (TT51) is valid also in the steady state, i.e. that 9 = 1 
and rj = 0. 
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We thus conclude that there exists a discontinuous 
transition at 7 = 0 + such that the Weibull distribution is 
swept in favor of a distribution that can support a finite 
flux at A = 0 . For completeness we discuss the evolution 
of the pdf upon straining in the Appendix. 

IV. CONCLUSIONS 

Our analysis of a generic 2 -dinrensional system indicated 
the existence of three types of solutions: 

(i) At 7 = 0 an isotropic Weibull solution 

p w ( A) = ^^(A/(A)) e exp-(A/(A)) 1 + e (22) 

exists which is singular as A —> 0 and for which there is 
no current; J = Oas0«2.2in the case of slow quench. 
In this state of mechanical equilibrium both the scaling 
law © and the scaling law m are obeyed. 

(ii) For any strain not equal to zero the statistics is fun¬ 
damentally different. For any value of the strain, both 
before and after the plastic yield there exists a finite flux 
of eigenvalues towards A = 0 where there is an absorbing 
boundary condition. This requires 0 = 1 . Both the scal¬ 
ing laws © and m appear to hold as long as 7 < y Y . 
Consequently we expect to find 77 = 0 in agreement with 
the numerical simulations. If the scaling law ([ 3 ]) is still 
obeyed, we expect to find /3 = — 1 . The direct numerical 
simulation found /3 ss — 0.92 (cf. Fig. [T| middle panel. 
We ascribe this and the smooth down fall to this value 
in Fig. [ 2 ] to finite size effects. 

(iii) Contrary to the analysis shown in the appendix 
of the breakdown of the Weibull distribution, we do not 
have a similar analysis of the yielding transition. It was, 
however, possible to solve for the steady state pdf in the 
form of Eq. (l 20 l) . This solution will be valid for 7 > 7 y- 
Indeed here 0=1 (compared to 0 = 2.2 in the isotropic 
case) and there is a steady state flux through the sys¬ 
tem. The actual shape of P ss (A) is not universal for large 
A and will depend on the form of the injection feeding 
jin (A). In Fig. [ 5 ] we have assumed constant feeding for 
all A < A max = 2 (A) . In contrast the shape for small A 
is expected to be universal. Since 0 = 1 we expect and 
find in this range 77 = 0 due to the ubiquitous scaling law 
m ■ On the other hand the scaling law © does not 
hold in the flowing steady state. This scaling law stems 
from the independence of plastic events for 7 < q Y and 
it fails in the steady state due to the build-up of correla¬ 
tions between the subextensive plastic events. We thus 
expect /3 = — 2/3 universally in this regime, in contra¬ 
diction with the scaling law ©. Our expectation for the 
theoretical dependence of /3 on 7 is thus shown in Fig. 0 

Finally, a word of caution is called for. As explained 
above in several occasions it is the distribution of min¬ 
imal eigenvalues of the Hessian that is the focus of our 
analysis. There is another distribution that may be at 
play here, i.e. the full distribution of plastic modes. This 



FIG. 6. Schematic presentation of the theoretical prediction 
for the 7 dependence of —(1 + 1 // 3 ). We reiterate that the 
value of this exponent at 7 = 0 is not universal, whereas in 
the steady state it is universal 

may have a different form and may be associated with 
other exponents. We assumed that the exponent /3 is as¬ 
sociated with the distribution of minimal eigenvalues. It 
remains to be seen whether one can determine the full 
distribution of plastic modes and how that distribution 
is related to the system-size dependence of (A7). At this 
point the full distribution of the eigenvalues of the plas¬ 
tic modes is not available for a generic amorphous solid. 
It will be certainly worthwhile to find new methods of 
determining it. 
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Appendix A: The Evolution of the Probability 
Distribution Below Yield 

In order to study qualitatively the evolution of the prob¬ 
ability distribution below yield we shall examine the sim¬ 
plest Fokker-Planck equation consistent with a singular 
drift velocity due to the saddle node dynamics under 
strain, namely 

dP/dj » ad[P/X]/d\ + D^)d 2 P/dX 2 +j in ( A; 7) (Al) 

and wish to study the evolution from the initial condition 
which is the isotropic Weibull distribution valid at 7 = 0 
and given by Eq. J 22 l Thus 

-P(A; 7 = 0 ) = P W (X). (A2) 

The easiest way to do this is to define the linear Liouvillc 
operator 


L = ad[l/X\/dX + Dd 2 /dX 2 


(A 3 ) 





7 


in terms of which we can rewrite Eq. (EH) as 

dP/dy = LP + j in { A; 7). (A 4 ) 

This equation has the general solution P = Phom + 
Pinhom where 

Phom( A; 7) = e yL P w ( A) 

PinhomiX-y) = e' 1,L [ e -7 L j in (X; y')dy'. (A 5 ) 

Jo 

The inhomogeneous solution depends on the structure 
of the injection rate ji n (A, 7) and is not universal. The 
only thing we really know is that the reinjection rate 
must ensure conservation of probability. We shall take 
the constant reinjection rate ji n ( A, 7) = ji n ( 7) between 
0 < A < A maa: in line with the avalanching process that 
occurs in the dynamics, while ji„(A, 7) = 0 otherwise, 
and consequently 

rl 

Pinhomi 7) = / *n(7W, (A6) 

Jo 

where ji n (y') is fixed by the requirement that probability 
is conserved. As P(A,7) = Ph om (A; 7) + P in hom{ y), we 
see that on integrating this expression between 0 < A < 
A max that 

r^max 

1= Ph om (A; y)dX + A max Pin horn (y)dX, (A 7 ) 

J y/2ary 

and this yields an expression for j in { 7) 


From Eq. ( 1 A 11 I) we see that as A —> yJ 2 ay from above we 
find 

Phom(y/^ay',y) = —j= / dtP w (X(y/ 4 Dyt) exp [-f 2 ] , 
v 7r JO 

( A1 2) 

because here A = 0 . But for A < ^/ 2 ay there is no solu¬ 
tion. Any nonzero value for P(A < 1/207,7) must come 
from feeding combined with flow from A = y/ 2 a 7 and 
adsorption at A = 0 resulting in the boundary condition 

P(A = 0 , 7 ) = 0 . 

On the other hand if A \f 2 wy, the Weibull term 
P w ( A(A, 7) + \Z 4 Dyt) can be expanded in powers of 
t inside the integral representation for Phom given by 
Eq. (IA 11 I) as the exponential exp — t 2 will ensure that 
only values of t <C 1 are important. In this manner an 
expansion in terms of the gradients of the Weibull distri¬ 
bution can be derived. To second order we find 

Pout (A; 7) = Oi (A; y)P w (A(A; 7)) 

+ C' 2 (A,7)Pu,(A(A; 7 ))' 

+ C 3 (A; 7 )Pui(A(A; 7 ))", (A13) 

where 

Ci (A; 7) = [1 + Erf( A(A; 7)/ y/ 4 D ^)]/2 
C 2 { A;7) = \fPhTn exp [-A(A; 7 ) 2 /( 4 P> 7 )] 

C 3 (A; 7 ) = (P7/2)[1-_ 

A(A; 7)/V (nDj) exp (—A(A; y) 2 / 4 Dy) 

+ Erf(\(X;'y)/y/ 4 D'y)], (A 14 ) 


-1 


jin( 7) = T- d{ 

^max 


' y/2aPy 


P horn (X]^)dX}/d , y. (A8) 


Note that the lower bound of the integral involving Phom 
only extends to ^ 2 ay. This (as we will see below) is 
because P? lom (A < 1/207; 7) = 0 . 

The homogeneous solution shows the evolution of the 
Weibull distribution under straining it and obeys 

/»00 

Phom (A; 7) = e yL P w ( A) = / dxP w {x)e lL 8 (x - A), 

Jo 

(All) 

The Liouville operator acting on the delta function 
e^ L S{x - A) = ^73= exp[~ (a ~^ 7)) ] will transform 
it into a Gaussian. In Eq. (IA 9 I) the expression A(A;7) = 
a/A 2 — 2 aq is the value that an eigenvalue A assumes due 
to drift alone. Thus we can rewrite Eq. ( 1 A 9 I) as 


Pftom(A, 7 ) — 


1 


\/4nDy J Q 


dxP w ( x ) exp 


(x — A(A; 7)) 2 


4 Dy 

(A10) 

Using the new variable t where x = A(A; 7) + \[ 4 .TTyt we 
can rewrite Eq. (lAlOll as 


P/iom(A,7) — 


\px 

y/AD-y 


dtP w (X{X\ y)+y/ 4 Dyt) exp [-t 2 ] 
(All) 


and the error function is given by Erf(x) = 
(2/1/7F) f Q exp [—t 2 \dt. Note once again that an examina¬ 
tion of Eo. lAldl shows that this solution becomes singular 
when A(A;y) = 0 or for A < \J 2 ay. In other words this 
solution is also only valid for A > y/ 2 ay. Thus for any 
7 / 0 a growing region exists where the Weibull solution 
breaks down. It is for this reason that we have written 
Phom = Pout in Eq. (IA 13 D to represent the outer solution 
of the distribution. 

We note there appears to exist a growing boundary layer 
with an inner solution P;„(A, 7) for A < y/ 2 ay. To solve 
for this inner solution we write 


P*n (A; 7) = E„A„( 7 )A" (A 15 ) 

and substitute in Eq. EH). For A < \J 2 ay, the Weibull 
distribution is swept away to to be replaced by essentially 
a constant flux distribution P,: n (A; 7) ~ Ai( 7)A for small 
A. In that case, by matching solutions at A = \J 2 ay we 
would find Ai{y) ss P out (y/ 2 ciy 1 7)/ \f 2 aJy due to fitting 
of inner and outer solutions. Together these arguments 
yield 


Pm (A; 7) w Pout {y/ 2 arf,y)/y/ 2 ay 


A. 


(A 16 ) 


As we can approximate Eq. (IA 12 I) by its saddle point 
approximation for small 7 we also find 


Phom ( l/ 2 ^ 7 ) 


( 29 Dy) 


[ e(xy- 


9/2 


4 nDy 


(A )(1 + (A »■ 


(A 17 ) 








































We are now in a position to calculate the flux J(7) and find that at small strain it grows in a singular manner as 

J( 7 ) ~ 7 0/2 - (A 18 ) 
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